arXiv:1503.00354vl [gr-qc] 1 Mar 2015 


Search for an emission line of a gravitational wave background 

Atsushi NishizaweEl 

Theoretical Astrophysics 350-17, California Institute of Technology, Pasadena, California 91125, USA 


Naoki Seto 

Department of Physics, Kyoto University, Kyoto 606-8502, Japan 
(Dated: March 3, 2015) 

In the light of the history of researches on electromagnetic wave spectrum, a sharp emission 
line of gravitational-wave background (GWB) would be an interesting observational target. Here 
we study an efficient method to detect a line GWB by correlating data of multiple ground-based 
detectors. We find that the width of frequency bin for coarse graining is a critical parameter, and 
the commonly-used value 0.25 Hz is far from optimal, decreasing the signal-to-noise ratio by up to 
a factor of seven. By reanalyzing the existing data with a smaller bin width, we might detect a 
precious line signal from the early universe. 


Newton refurbished our understanding of the colors of 
visible light by sophisticatedly utilizing prisms, in addi¬ 
tion to his extraordinary accomplishments on gravitation 
and classical mechanics. His sequential works on optics 
were summarized in Opticks published initially in 1704. 
But we had to wait for another ~ 100 years for signifi¬ 
cant experimental developments on electromagnetic wave 
spectra. In 1800, Herschel uncovered infrared compo¬ 
nents of a light ray. Around the same time, Wollaston 
(more systematically Fraunhofer) discovered sharp ab¬ 
sorption lines in the spectrum of Sun light. Moreover, 
the new element Helium was identified in the solar corona 
by its emission lines detected at the total solar eclipse in 
1868. However, the physical mechanism itself behind the 
observed sharp emission/absorption lines could be hardly 
clarified in 1800s. Indeed, this conundrum was one of the 
primary driving forces toward the construction of quan¬ 
tum mechanics in early 1900s. 

On the other hand, the existence of a gravitational 
wave (GW) is predicted by Einstein’s general theory of 
relativity that was completed in 1915. As the gravita¬ 
tional interaction is very weak, GWs have not been di¬ 
rectly observed yet. However, with the advent of power¬ 
ful second generation interferometric detectors, the first 
detection would be achieved shortly. 

The high penetrating power of GWs could become 
quite advantageous to probe the early universe. At 
present, our observational information before the big 
bang nucleosynthesis is severely limited [l|. Therefore, 
once primordial GWs are detected, they would provide 
us with invaluable clues to understand inflation and the 
subsequent reheating era [2-0 and high energy physics 
0H3- Also GW backgrounds (GWBs) bring us oppor¬ 
tunities to test gravity theories [1 ll4l4|. To maximally 
extract the scientific outputs from the accumulated data 
of the interferometers, we should deliberate on meth¬ 
ods of data analysis and thoroughly search for stochastic 
GWBs, not only for given theoretical predictions but also 
in model independent manners. 

In the light of the aforementioned history of electro¬ 


magnetic waves, it would be interesting to explore a sharp 
emission line of a cosmological GWB. In this letter, we 
study how a line GWB looks like in detector signals, and 
suggest an efficient detection method by correlating data 
of multiple ground-based interferometers. We find that 
the width of frequency bin should be chosen carefully, and 
its preferred size is smaller than 0.1 Hz. With the com¬ 
monly used width 0.25 Hz [lj|, the signal-to-noise ratio 
could be reduced by up to a factor of seven. Therefore, by 
reanalyzing the existing data, we might strike a precious 
line signal that has been buried in detector noises. 

Doppler broadening of GWs — Throughout this letter, 
GWB is assumed to be isotropic in the CMB rest frame. 
As a simple but suggestive example, we first examine 
a delta-function-like spectrum concentrating solely at a 
frequency / r (cosmologically generated at the same co¬ 
moving wavelength). Our study can be easily extended 
to a GWB with multiple lines. The effects of a line profile 
will be discussed later. 


For a detector moving relative to the CMB frame with 
velocity v(t), the observed Doppler-shifted frequency 


f(,t,Q) = / r 1 - v{t) ■ n 


is anisotropic and depends on 


the propagation direction of a GW O. In this letter, we 
use the unit c = 1 for the speed of light, and also drop 
the terms of 0(|f7(t)| 2 ) for the Doppler effect. 


The relative velocity v s of the Solar-System barycen- 
ter (SSB) to the CMB rest frame is toward the Galactic 
coordinate {£,b) = (263.99° ± 0.14°, 48.26° ± 0.03°) with 
the magnitude v s = |F S | = (1.230 ± 0.003) x 10 -3 (corre¬ 
sponding to 369.1kms -1 ) [16j- Below, for simplicity, we 
put v(t) = v B , neglecting corrections (at most <10%) due 
to the velocity of the detector relative to the SSB. We also 
neglect minor relativistic effects such as the Sachs-Wolfe 
fluctuations. These would be reasonable approximation 
for our demonstration, and we can easily employ more 
sophisticated models in actual data analysis. Then we 
have 


f(t, fi) = /r [1 — v 8 cosd] 


(1) 
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FIG. 1: Observation of a line GWB that is isotropic in 
the CMB rest frame. An observer (an interferometer on 
the Earth) moves at the relative velocity v(t), and sees an 
isotropic background due to the Doppler effect. The ob¬ 
served frequency is identical along the ring across the sky 
with 9 = const (8: angle between v(t) and the propagation 
direction of GW fl). The spin rotation of the Earth generates 
diurnal modulation for the correlated signal of two interfer¬ 
ometers I and J. 


with the angle 6 between 12 and v s . For a given observa¬ 
tional frequency, the GW signals come from a ring on the 
sky as shown in Fig. |T| The total width of the Doppler 
broadening is given by 

V D=2/.„ s ~0.2( Ti A i -)H, (2) 

The optimal band of ground-based interferometers is 
around 10-1000 Hz. Meanwhile, by taking Fourier trans¬ 
formation for an appropriate time-segment (duration T ) 
of data, we can realize a fine frequency resolution ~ 
10 _3 (T/1000 sec) -1 Hz, unlike an electromagnetic-wave 
spectrum. But, as we discuss later, the spin rotation of 
the Earth plays an important role and this limits the 
duration T for the coherent Fourier transformation. 

GW signal — Next, we briefly summarize the correla¬ 
tion analysis for an anisotropic GWB (not specific to a 
line GWB) E3- In the frequency space, the GW signal of 
a detector at the position Xj can be generally expanded 
as 

Hf , */) = V f dCl h A (f , 0) e - 2 -^'^ Ff{Cl) . (3) 

Here, the label A denotes the polarization mode (A = 
+, x) of a GW, and h A {f, Cl) is the Fourier mode of the 
GW in each polarization state. The factor Ff is the an¬ 
tenna pattern function of the J-th GW detector, and has 
a quadrupolar angular dependence. The power spectrum 
Sh(f) for h A (f, Cl) is given by the energy density of GWB 


H gw (/; Cl) (per logarithmic frequency interval per stera- 
dian) normalized by the critical density of the universe 
as 


A 47T^ 

Cl sw (f-,Cl) = ^fS h (f;Cl), (4) 

where Hq is the Hubble parameter and we use Hq = 
70 kms -1 Mpc -1 . 

The correlation analysis is an efficient method to dis¬ 
tinguish a GWB signal from detector noises (see e.g. 0 - 
13). In this method, we prepare data streams of two 
widely-separated detectors (with independent noises), 
and take cross correlations of their Fourier modes. For 
statistical amplification and coarse graining, we divide 
the Fourier space into multiple bins. Below, a bin is spec¬ 
ified by the label k , and its central frequency and width 
are denoted by fk and Sfb , respectively. We intrinsically 
have Sfb > T -1 due to the frequency resolution. Within 
each bin k , we take a summation of the correlated prod¬ 
ucts. Its expectation value is written as [20j 

~ 87r 2 / 3 Z ^ k ^' 

Here we defined the integral 


Z{h) = 


dQ 


^gw {fk\Cl)l{fk,Cl) 


( 6 ) 


j 4tt gw 

and the direction-dependent overlap-reduction function 


lift ^) = ex P 27r */^ • du 


( 7 ) 


with djj(t ) = Xj{t)—Xj(t) for two detectors I and J. At 
high frequency regime, the phase factor exp ■ du 


depends strongly on the propagation direction Cl. In 
Eq. ©, we assumed to take a small bin width Sfb so that 
the frequency dependence of relevant functions can be 
neglected within each bin. But this is not always the 
case, as we later discuss in detail. 

Now we apply the above expressions specifically to the 
line GWB at a frequency / r (defined in the CMB frame). 
From Eq. ©, we can write the observed spectrum by 


^gw(/; d ) = g gw/r $[f ^ /r(l - v s cos (9)], (8) 

where e gw is the total energy of the line normalized by 
the critical density. 

For a line GWB, we define the over lap-reduction func¬ 
tion F(/, u) that is obtained by integrating y(/, H) across 
the ring directions at u = cos 9 = constant as 

r (/>) = ^ ( 9 ) 
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This definition is different from the following standard 
one 7 s t(/) by an additional u- integral and the overall 
normalization; 

r rl p2tt 

7st (f) = —j du J #7(/,0). (10) 

From Eqs. ©-© we have 

Z{h) = ^ J 1 du S[f k -f r { 1 - « B u)] r(/ fc) u) . 


Replacing the frequencies of slowly-varying functions 
with / r and putting T(f k ) = F(/ r , «(/*,)), we obtain 


_ ^e gw J( f ^ 

Mfc — 2 r3 ' 

lb7T Z 


( 11 ) 


Note that the frequency / and the sky angle 9 (equiva¬ 
lently u) are related by Eq. ©■ 

In Fig. [21 we show the directional overlap-reduction 
functions T(fk) observed by the two advanced-LIGO 
(aLIGO) interferometers for a line GWB at / r = 200 Hz. 
The frequency dependence changes diurnally due to the 
spin rotation of the Earth. 

The left and right edges of the curves are for GWs 
coming from the cold and hot spots of the CMB dipole, 
and have large amplitudes. This is because the GW sig¬ 
nals are hardly canceled by the 0-integral in Eq.©. In 
Fig.2, we should also notice that the real (imaginary) 
part is symmetric (anti-symmetric) around the intrinsic 
frequency / r . In relation to this structure, for an isotropic 
GWB ( e.g. without a line component), the imaginary 
parts of the correlated signals vanish due to a cancella¬ 
tion. In contrast, for a line GWB, the overlap-reduction 
function T(/) can have finite imaginary parts because of 
the anisotropy induced by the Doppler effect. 

In Fig.2, the wavy structure of the function T(f) is 
mainly determined by the phase factor in Eq. ©. An 
oscillation occurs when the phase changes by 8a ~ 2-7T, 
where a = 2nf r djj cos /3 and cos (3 = D • dij. We also 
have <5(cos/3) ~ <5(cos0) ~ l/(/ r d/j). Then the charac¬ 
teristic frequency interval Sf c for the oscillation is evalu¬ 
ated as 

~ /.o,«(cosffl ~ - 0-1 '=-■ (12) 

Comparing this result with the Doppler width ©, the 
number of oscillations becomes larger for a higher fre¬ 
quency / r . For another detector pair more separated than 
aLIGOs (du ~3000km), the characteristic interval 8f c 
becomes smaller. 

Signal-to-noise ratio — The correlated signal m 
would be suitable for a matched filtering analysis with 
the two fitting parameters / r and e gw . To calculate the 
signal-to-noise ratio (SNR) of the line detection, we eval¬ 
uate the variance of the correlated signal n k . Here, for 



/ 


FIG. 2: Overlap-reduction function F(/) = F(/ r , u (/)) for two 
LIGO interferometers with / r = 200 Hz. The observed fre¬ 
quency / and the propagation direction u = cos 6 are related 
by Eq. Q. The upper panel is the real part and the lower 
panel is the imaginary part. The overlap-reduction function 
is modulated by the spin rotation of the Earth. Its shape is 
plotted at four different epochs; the fiducial time (blue, solid), 
3h later (orange, dashed), 6h later (green, dotted), and 9h 
later (red, dot-dashed). The characteristic frequency interval 
Sfc ~ 0.1 Hz of the oscillations is given by Eq. (Ell- 


each Fourier mode, the noises of the two detectors are as¬ 
sumed to be independent and to have magnitudes much 
larger than those of the GW signals. Then, the variance 
of the correlation signal is mainly contributed by de¬ 
tector noises and is given by [2pj 

<d*j6f b P I (f k )Pj(f k ) , (13) 

with the one-sided noise spectra Pi(f) and Pj{f). 

The correlation signal in Eq. © is a complex number, 
as mentioned before. From Eqs. ©, CD, and m, the 
total squared SNR is evaluated as 

2 l^fcl 2 _ 9 -Hq Te^Sfb 2 

P 2-, a 2 ~ 32-71-4 f®v% Pl (f r ) Pj (f r ) ' 

(14) 

For a total observational duration T 0 b s much longer than 
a day, we need to divide the data into time segments for 
taking Fourier transformation. But, in Ea. (114L we can 
effectively regard T as the duration T 0 b s , after replacing 
|F(/fc)| 2 by its time average. Importantly, the total SNR 
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p does not depend on the details of the time segmen¬ 
tation, once we can take a sufficiently small bin Sfb to 
resolve the function r (/*.). 

Signal cancellations within bins — So far, we have as¬ 
sumed that the bin width Sfb is much smaller than the 
characteristic frequency interval Sf c . If this is not the 
case, the coarse-grained signal \pk\ decreases due to the 
averaging of the wavy structure around zero point, as 
understood from Fig.2. Then, for Ea. (ll4l) . we should use 
the following integral 

3Hq Te gw f jj-tv f \ 

Mfc = i ~ 2 . 3 / dfT(f), (15) 

16tt 2 ffv s J 5fb 


in stead of the original one m that is valid only when 
the variation of the function T(f) is negligible within each 
bin. Therefore, to realize a high SNR, we should set 
the bin width sufficiently smaller than the characteristic 
interval Sf c ~ 0.1 Hz. Nevertheless, the commonly used 
value is Sfb = 0.25 Hz, e.g. in [l5|, and the sensitivity to 
a line GWB might be decreased significantly. Below, we 
take a close look at this degradation. 

Sensitivity to a line GWB — In Fig. [31 the sensitivity 
to the normalized energy density e gw is shown as a func¬ 
tion of the line frequency / r . Here, we assumed a lyr 
observation with two aLIGO detectors at the detection 
threshold of p = 10. For their noise spectra P/,j(/), we 
use the fitting formulas given in [21]. 

In our study, the key parameter is the width of the fre¬ 
quency bin Sfb . For our demonstration here, we selected 
the following three widths; 0 Hz, 0.25 Hz, and 2.5 Hz. The 
first one means a sufficiently small width (still satisfying 
Sfb > T~ l ). We use Eq. (fT4l) for this one and Eq. (fL5l) for 
the other two. Because, for a finite width, the sensitivity 
depends on the actual positions of the bins, we plotted 
the sensitivity averaged over their possible positions (by 
shifting boundaries of bins). 

In Fig.3, the sensitivity decreases significantly for 
larger bin widths. This can be expected from the previ¬ 
ous arguments on the cancellation of the wavy structure. 
The degradation of the sensitivity, say, at / r = 70 Hz, is 
by factors of 6.8 for Sfb = 0.25 Hz and 22 for 2.5 Hz. Note 
that these ratios are independent of the noise spectra. We 
also examined bin widths smaller than Sf c =0.1 Hz. For 
a line GWB with 10Hz < / r < 600Hz, the degradations 
(compared with Sfb = 0) are within 1.1 for Sfb = 0.01 Hz, 
1.2 for 0.02Hz, 1.5 for 0.04Hz, and 3.5 for 0.09Hz. For 
example, to keep the loss of the sensitivity within 20%, 
we should take a bin width Sfb smaller than 0.02 Hz, in 
contrast to the standard choice 0.25 Hz. 

In the opposite limit with Sfb <5/d, the /-integral in 
Eci. ([T5l) is equivalent to posing an additional u-integral 
to Eq. ©• The resultant is proportional to the standard 
overlap-reduction function 7 s t(/) defined in Ea. fXOl) . In¬ 
deed, in Fig.3, the positions of the peaks are around the 
zero points of the standard one 7st(/) (20j . 



FIG. 3: Sensitivity to the normalized energy density e gw of a 
line GWB. We assumed a 1 yr observation with two aLIGO 
detectors and the detection threshold of p = 10. The lower 
curve (blue, solid) is for the bin width 5b = 0 Hz. The middle 
(orange, dashed) and upper (green, dotted) are for the bin 
widths of 0.25 Hz and 2.5 Hz, respectively. 


Discussions — In the analysis above, we assumed a 
delta-function-like profile for a line GWB. Here we briefly 
comment on the effects of the intrinsic line width A/ gw 
(now only considering Sfb <C Sfc)- This width A/ gw 
should be compared with the Doppler broadening Sfr> 
and the characteristic frequency interval 6f c . For a sharp 
line with A/ gw <C Sf c , we can apply our previous argu¬ 
ments without change. 

In the limit of a broad spectrum with A/ gw <5 /d, a 
GWB at a certain observed frequency is contributed from 
the whole sky directions by the Doppler effect. There is 
no correspondence between an observed frequency / and 
the angle 9. In this case, by replacing H gw (/fc;H) with 
figw(/fc) and using Eqs. ©, (flOl) and (TUH) . we have 


_ 9^0 rp r r ^gw(/fc)7st(/fc) 

"50^ ^ f[]Pi(fk)Pj{fk) ' 


(16) 


This is the standard SNR formula for a broad-band 
isotropic GWB J20[. Note that, for a fixed total en¬ 
ergy density / dln/H gw (/) = constant (e.g. for the 
bound from the big bang nucleosynthesis [l] or CMB 
EH), we can readily derive an approximate scaling re¬ 
lation p ex (A/gw) -1 / 2 for the SNR. 

In intermediate case with Sf c < A/ gw < Sfo, a GW 
signal at a certain observed frequency is contributed by 
GWs from a band region in the sky. Thus the GW signals 
are partially canceled, and the SNR is smaller than that 
for A/ gw = 0 (for a fixed energy density). 

When searching for a line GWB with real detectors, 
one might worry about confusion with various noise lines 
produced by instruments. However, we can clearly distin¬ 
guish them by using their time dependences. In contrast 
to noise lines, a GW line is modulated due to the rotation 
and revolution of the Earth, showing the characteristic 
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profiles as illustrated in Fig.2. Therefore, by appropriate 
Matched filtering analyses, one could easily distinguish 
the true GW signals from artificial noises. 

Although we have only considered aLIGO pair, our 
study can be also applied to other planned detectors, such 
as LISA (sensitive around 10 -4 -10 -1 Hz) and Einstein 
Telescope (ET) (1-1000Hz). With the original LISA of 
the three-arms configuration [23], we can construct two 
noise-independent data streams [ 3 ]. But, due to the 
symmetry of the system, their correlation identically van¬ 
ishes for an isotropic GWB (as usually assumed). In 
contrast, a line GWB appears as a strongly anisotropic 
signal and could be detected, in principle, by the corre¬ 
lation. The situations are similar for ET with the pro¬ 
posed three-arm configuration jHI- However, eLISA is 
currently planned to have just two arms and can generate 
only single data stream [26|. In this case, our correlation 
method no longer works. For a space mission like LISA, 
its third arm is essential to explore a line GWB. 

Conclusions — We have studied the method for effi¬ 
ciently searching for a line GWB. The sensitivity might 
be degraded significantly, unless the bin width 5ft, is 
much smaller than 0.1 Hz (in the case of two LIGO pair). 
Nevertheless, for correlation analysis of GWBs, the com¬ 
monly used width is 0.25 Hz. By reanalyzing the existing 
data with a smaller width, we might actually uncover an 
important signal from the early universe. 
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